EXCEPTIONAL CIRCLES OF RADIAL POTENTIALS 
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Abstract. A nonlinear scattering transform is studied for the two-dimensional 
Schrodinger equation at zero energy with a radial potential. First explicit 
examples are presented, both theoretically and computationally, of potentials 
with nontrivial singularities in the scattering transform. The singularities arise 
from non-uniqueness of the complex geometric optics solutions that define the 
scattering transform. The values of the complex spectral parameter at which 
the singularities appear are called exceptional points. The singularity for- 
mation is closely related to the fact that potentials of conductivity type are 
"critical" in the sense of Murata. 
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1. Introduction 

We study singularities of the scattering transform at zero energy for two-dimensional 
Schrodinger operators with radial and compactly supported potentials. We will 
present what is, to our knowledge, the first example of a family of potentials for 
which the singularities of the scattering transform can be computed explicitly. The 
singularities, occurring at so-called exceptional points of the potential, arise from 
non-uniqueness of complex geometric optics (CGO) solutions of the Schrodinger 
equation. We also present the numerical computations which inspired this work 
and display the formation of singularities under perturbation of a potential of con- 
ductivity type. 

To motivate our main results, let us define a class of potentials that plays a 
central role in our work. 

Definition 1.1. A compactly supported real-valued potential q € C^°(IR 2 ) is of 
conductivity type if q — tp (A?/j) for some real-valued ip 6 C' X) (M. 2 ) satisfying 
ip(z) > c > for all z in a bounded set Q C K 2 and tp(z) = 1 for all z£K 2 \!l. 

Note that the positive function tp above solves (—A + q) ip = 0. This terminology 
arose when Schrodinger scattering theory was used to analyze the inverse conduc- 
tivity problem in Nachman [18] , and was also needed in [T4j [19] . In those works q is 
not necessarily compactly supported, but a condition implying hrri.| z |_>. 0o ip(z) = 1 
is crucial. In Appendix [C] we show that the associated Schrodinger operator has 
no eigenvalues, and that the function tp representing q is unique. 

Recall that the inverse conductivity problem of Calderon [5] consists in recon- 
structing the conductivity a of a conducting body C M 2 from the Dirichlet to 
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Neumann map, denned as follows. Let / £ ^^(dVl) and let u £ solve 

V ■ (o-Vti) = 
u \on = /■ 

This solution is unique if a £ L°°(Cl) is real- valued and strictly positive. The 
Dirichlct-to-Ncumann map is the map A CT : iJ 1 / 2 (9r2) — > H~ x / 2 (d£l) given by 

Kf = <r — . 
9l/ an 

Calderon's problem is to reconstruct the function a from knowledge of the map A CT . 
In dimension two, Calderon's problem in its original form was solved in [3]. 

Nachman [T5] exploited the fact that ip = a x l 2 u solves the Schrodinger equation 
(—A + q) ip = where q = o~ x l 2 A (c 1 / 2 ) . The Schrodinger problem (—A + q) tp = 
also has a Dirichlet-to-Neumann map A g which can be determined from A a . Using 
A g , Nachman was able to construct CGO solutions to the Schrodinger equation, 
use these solutions to compute the scattering transform of q, and reconstruct a by 
inverting the scattering transform. 

The theory of the scattering transform has been worked out in detail for po- 
tentials of conductivity type (see Lassas, Mueller and Siltanen [13] and references 
therein), but far less is known about general classes. Grinevich and Novikov [7J 
part I of supplement 1] commented that Schrodinger potentials with well-behaved 
scattering transforms are "not in general position," and Nachman proved under 
minimal hypotheses that the scattering transform is regular if and only if the po- 
tential is of conductivity type. Below we show how work of Murata |17j implies 
that the set of conductivity type potentials is unstable under perturbations and 
is therefore unstable in any reasonable function space! 

Let us define the CGO solutions and scattering transform for a potential q £ 
Cfi°(R 2 ). To define the CGO solutions, let k £ C, set z = x + iy and write 

(1.1) kz = (ki + ik?) [x + iy) 
(complex multiplication). The CGO solutions ip(z,k) solve 

(1.2) (-A + g)^ = 0, 

lim e- lkz ^{z,k) = 1. 
\z\— >oo 



The set of nonzero k £ C for which (1.2) does not have a unique solution is called 



the exceptional set £. If the exceptional set is empty, the solutions ip(z,k) form a 
smooth family, and the scattering transform of q is given by 

(1.3) t(k) = J e^q(z)ip(z,k) dz. 

The behavior of t at k = plays a special role which we will elucidate in what 
follows. 

Nachman showed (under rather less stringent regularity assumptions than q £ 
C^°(R 2 )) that q is of conductivity type if and only if: 

(i) £ is empty, and 

(ii) |t(fc)| < C|fc| e for \k\ small and some e > 0. 

(see [T5], Theorem 3). Until this time it was not clear how the scattering transform 
behaved for potentials outside this limited class. Our purpose here is to construct 
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and analyze examples for which the singularities of t(k) can be computed explicitly. 
In our examples, t(fe) is well-defined except on a circle in the complex fc-plane, and 
we can compute the singularity explicitly. 

We will study the scattering transform for families of radial potentials q\ defined 
as follows. Denote by B\ the open unit disc in R 2 centered at 0, so that dB\ = S 1 
regarded as an embedded manifold in R 2 . Suppose that a g C°°(]R 2 ) is a real- 
valued radial function satisfying a(z) > for all z <E K 2 and a — 1 <E (-Bi). 
Then go : = cr _1 / 2 Acr 1 / 2 e (Bi) is a radial potential of conductivity type. For a 
nonnegative radial function w € C^°(Bi) not identically zero, set 

(1.4) q x = q + Xw. 

We show in Appendix [A| that , for a radial potential, t(fc) is a real- valued and radial 
function. Our main result is: 

Theorem 1.2. Denote by t\ the scattering transform of q\. For small A =/= 0, 

log 

as k — > 0. Moreover: 

(1) For A > sufficiently small, the exceptional set £ is empty, and the scattering 
transform t\ is C°° away from k = 0. 

(2) For A < sufficiently small and a unique r(A) > 0, the exceptional set £ is a cir- 
cle C\ of radius r(A) about the origin, and the function t\ is C°° on M. 2 \ [C\ U {0}], 
while 

lim |tx(fc)| = oo. 
IfcKr(A) 1 

The radius r(A) obeys the formula 

(1 + (A)) 



2tt^(A) 



(1.5) r(A) ~ q cxp 2vr I ft 

w/iere 

2tt' 

7 is Euler's constant, and /x(A) is i/ie eigenvalue of the Dirichlet-to-Neumann op- 
erator for qx corresponding to the constant functions on S 1 . 



Theorem |1.2| shows that A = is an 'essential singularity' for the map A ^ t^. 
The cases A = 0, A < 0, and A > may be characterized in the following way. We 
recall from Murata {TT that a Schrodinger operator — A + q is called 

(i) subcritical if —A + q has a positive Green's function, 

(ii) critical if —A + q does not have a positive Green's function, but the qua- 
dratic form 



(v,v)= f (\{Vv){z)\ 2 + q{z)Hz)\ 2 ) 



dz 



on C^M 2 ) x ^(R 2 ) is nonnegative, and 
(iii) supercritical if the quadratic form q is not nonnegative. 

Appendix[B]below shows how [TTj implies the following. First of all, the conductivity- 



type potential qo is critical. Furthermore, taking A < in ( 1.4 1 gives a supercritical 



potential q\ which cannot be of conductivity type since there is no positive solution 



of (— A + q) ip = 0. Finally, taking A > in (1.4) gives a subcritical potential q\ 



allowing a unique positive solution tp of (—A + q) ij) = 0. However, this ip(z) grows 
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logarithmically in \z\ and does not satisfy limui^oo tp(z) — 1, so q\ is not of conduc- 
tivity type in the sense of Definition |1.1| Therefore, the class of conductivity-type 
potentials is not stable under perturbations of A. 

The scattering transform T '■ q — > t plays an important role not only in Nach- 
man's solution of the inverse conductivity problem, but also in the solution of 
the Novikov-Veselov (NV) equation by the method of inverse scattering (see Las- 
sas, Mueller, Siltanen and Stahel [H] and Perry [TH] for details and further ref- 
erences). Thus, a clear understanding of the singularities of T is important both 
for Schrodinger inverse scattering and for understanding the dynamics of the NV 
equation and other completely integrable equations in the NV hierarchy. 

It follows from Perry [19] that the Novikov-Veselov equation with initial data go 
has a global solution. On the other hand, Taimanov and Tsarev [24, 25, 26, 27J have 
used the Moutard transformation to construct initial data q for the Novikov-Veselov 
so that the solutions blow up in finite time. These potentials have L 2 eigenvalues 
at zero energy, and hence are not of conductivity type (see Proposition C.ll. Wc 
conjecture the following dichotomy for the Cauchy problem for the Novikov-Veselov 
equation: 

(a) If the initial data is critical or slightly subcritical, the solution exists globally 
in time, and 

(b) In all other cases, the solution blows up in finite time. 
We hope to return to this question in a subsequent paper. 

We close this introduction by sketching the proof of Theorem |1.2| and summa- 
rizing the contents of this paper. To analyze the singularities of the scattering 
transform for the family (1.4 1, we recall the reduction of (i) the problem (1.2) 
and (ii) the map 
type (see 



|1.3| ) respectively to (i) a boundary integral equation of Fredholm 
and (ii) a boundary integral (see 1.12). We refer the reader to 



Nachman 



Theorem 5 and its proof in section 7 for a complete discussion. 
In order to state the reduction, we first define the Dirichlet-to-Neumann map, 
denoted A q , for the Dirichlet problem 

(1.6) (-A + q)u = 0mB 1 

u lsi = /• 

If zero is not an eigenvalue of th e operator — A + q with Dirichlet boundary con- 
ditions on Bi, the problem (1.6) has a unique solution u for given / G H 1 ^ 2 (S 1 ), 
and we set 

(1.7) A q f °" 



dp 



where d/dv denotes differentiation with respect to the outward normal on S 1 . We 
let (•, A q -) denote the corresponding bilinear form: 



(1.8) 



<<?,A 9 /> 



(Vu • Vm + qvu)dz 



D 



v~dS 
ov 



We will denote by A 
It is known that A„ : 



where u solves (1.6) and v 6 H 1 (D) with v\ gl = g. 
the Dirichlet-to-Neumann operator for (1.6) with q = 
H 1 / 2 (dn) — > H^ 1 / 2 (dVl). If ip denotes the restriction of the unique solution of 
Ol to S\ then 



(1.9) 



Tip 
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where T : H 1 / 2 (S 1 ) -> H 1/2 (S r ) is the compact operator 
(1.10) Tifj = Sk (A g — A ) ip. 

Here, the operator Sk is an integral operator 



(1.11) 



{S k i,) (*) = / G k (z - z')^{z') dS(z'), 
Js 1 



Gk( ■ ) is Faddeev Green's function (see (2.4) in what follows), and dS is arc length 
measure on S 1 . The formula 



(1.12) 



t(fc) = / 
J s 1 



, / k z 



[(A,-A )^] (z,k) dS(z) 



computes the scattering transform in terms of the boundary data tp that solve ( 1.9 1 



To prove Theorem 11.21 we will study the family of operators T^a corresponding 



to (1.10) with q given by (1.4). We will show that the resolvent (J + T^ ^) has 
a rank-one singularity if A < for |fc| = r(A), where the asymptotics of r (A) are 



given by ( 1.5 ). We will then use ( 1.12 ) to show that this rank-one singularity leads 
to a singularity in t (k) on the circle of radius r(X). 

The structure of this paper is as follows. In jj2rwe recall some basic facts needed to 
analyze ( 1.9 1 and the associated operators. In q3]we extract the rank-one singularity 



of the resolvent (I + T)^ 1 and show the existence of a 'circle of singularities' for t\. 
Finally, in f|4j we present numerical computations of and compare the results to 
the analysis in |3j 
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2. Preliminaries 

First, we recall that i? 1 / 2 ^ 1 ) can be defined in terms of the Fourier basis {fn} n —- 
of ^(S 1 ) given by 



(2.1) 



<Pn(0) 



/2tt 



in the following way. A function / g L 2 {S X ) belongs to i? 1//2 (S' 1 ) if / = J2 n a nfn 
and 



(2.2) 



1/11^1/2(51) = ( X + H) I ™'" 



n— — 00 



is finite. If we denote by P the projection 
(2.3) ■ ' 



(PVO (z) = — ij{w) dS(w), 

27T Jgl 
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and let Q = I — P, then P and Q are orthogonal projections in 7J 1//2 (5 1 ). We 
denote by H^ 1 / 2 {S 1 ) the topological dual of iJ 1 / 2 (S' 1 ) and by ( ■ , • ) the inner 



product associated to the norm (2.2) 



Next, we recall some basic facts about the operators Sk and A q that appear in 



the boundary integral equation (1.9) 



First of all, we need some properties of the integral operator Sk ■ A reference for 
this material is the thesis of Siltanen [3T]. For k € C, the operator Sk in ( |1.11[ ) is 
defined in terms of Faddeev's Green's function (setting z = x + iy) 

(2.4) Gk(z) = e lkz g k (z) 



where kz is given by ( 1.1 1 and 



(2-5) 9k (z) = 7^/ t ~ ''^ dtidb 

(2tt) 2 J r2 £ (£ + 2k) 



Here £ • z — £ix + In the denominator of ( |2.5[ ), £ = £1 + i&- The function gk is 
the convolution kernel for Green's function for the operator —i (<9 (d + fc)J . The 
factor e lkz normalizes G k to be a fundamental solution for A = Add. Thus 

(2.6) H k (z) = G k (z) - G Q (z) 
is smooth and harmonic, where 

G (z) = ~]og\z\ 

is the normalized fundamental solution for —A. The integral operator Sk is a 
bounded operator from iJ _1 / 2 (S' 1 ) to il 1 / 2 (5 1 ) (see, for example, Lemma 7.1 of 
Nachman [18), but we will need a finer description. The following Lemma is a 
simple consequence of [21], Theorem 3.2 and following discussion), and we omit the 
proof. 

Lemma 2.1. The formula 

(2.7) Hk(z) = H 1 (kz)~^-\og\k\, 
holds. Here H±( ■ ) is real-valued, smooth, and harmonic, and 

where 7 is Euler's constant. 

From the Lemma, we immediately conclude: 
Lemma 2.2. The decomposition 

(2.8) S k = S Q +U k -(\og\k\)P 
holds, where 

(2.9) («fcV0(*)=/ ffi^-^VsfA 

Remark 2.3. By straightforward computation, SqP = so that 

(2.10) (SfcPVO (z) = [2*Hi{kz) - log |*|] (PV>) («) 

where we have used the mean value property for harmonic functions, while 

(2.11) SkQ = S Q Q + nkQ. 
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We will use this decomposition to analyze the singularities of the operator Tk,\ as 
A tO. 

Next, we need some simple properties of the Dirichlet-to-Neumann map and the 
integral operator T in the presence of radial symmetry. First, if q is radial, the 
problem (1.6) can be solved by Fourier analysis on the circle using the basis (2.1). 
We set, for z — re lB , 

oo 
n— — oo 

To compute the Dirichlet-to-Neumann map, we solve the problem 

(-A + q)i> =0 in B x 

i>\ S i =^Pn 

Writing ip = ^ n (r)ip n {9) we have 

r dr \ dr J ' \r 2 

It follows that 

(2.13) A q (p„ = fj, n (q)fn 



•Ml) 



where 



Thus A g has a complete set of orthonormal eigenf unctions and real eigenvalues, and 
hence commutes with complex conjugation. Furthermore, it is easy to see that 

(2.14) = n n G K. 
Let 

(2.15) M (A) := /xo(«a), 



where q\ is given by (1.4). The following fact is crucial. 

Lemma 2.4. Suppose that a G C°°(IR 2 ) is a real-valued radial function satisfying 
a(z) >0forallz£ M 2 and a - 1 G Cg° (Si). Denote q := a^Aa 1 / 2 G Cg° (Si). 
Let u) G C^°(-Bi) fee a nonzero, nonnegative radial function, and define q\ by (1-4-). 
Then fi(0) = and //(0) > 0. 

Proof. The function cr 1 / 2 is the unique solution to (—A + go) 4> — with <j>\ sl = 1, 
and since a 1 / 2 is constant in a neighborhood of the boundary we have -§^ (jl ^ 2 \ s i = 
0. It follows that /Lt(0) = 0. The fact that /x(A) is continuously differentiable in A 
follows from the fact that the unique solution to the problem (2.12) with q given 
by (1.4) depends analytically on A. 

To compute /i'(0), let ip\ solve the Dirichlet problem 

(2.16) (-A + <z A )^ A = 0, 

V'Algi = 1, 
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and let ip\ = dijj\/d\. Since /i(A) = (1,A 9A 1), we calculate fi'(0) by taking the 
derivative of this pairing: 

/,'(A) = ^<1,A^1> = ~ / \V^\ 2 +q x \^\ 2 dz 



ID 

" v W\ T ^q\WW\ T w\H>\ ' " ' 

id 

= 2(^ a | s i,A 9a 1)+ / w\^ x \ 2 dz. 

Jd 



We have ip> 



^ = 0, so //(O) = J D w\ip \ 2 dz > 0. 



□ 



3. Singularities of CGO Solutions and Scattering Transform 



In this section we compute the singularities of solutions to the boundary integral 
equation ( |1.9| . Our computation is based on the following decomposition of the 
resolvent (I + T ki \) 1 as a bounded operator from H 1 / 2 {S 1 ) to itself. We note 
that, since qg is of conductivity type, it follows from Theorem 5 in |18j that the 
operator (J + T k0 ) has a bounded inverse for all k 6 C. 

Lemma 3.1. Suppose that a € C°°(IR 2 ) is a real-valued radial function satisfying 
a{z) > for allze M 2 and a -I e C§° {B x ). Denote q := a^^Aa 1 / 2 G Cg° (Bi). 



Let w £ C^°(Bi) be a nonnegative and nonzero function, and define q\ by (1-4-) 



Suppose that T k \ is given by (1.10) with q = q\. Then, for all |A| sufficiently small, 
there is a rank-one operator F = F(k, A) and a bounded operator R = R(k, A) so 
that 

(3.1) (J + Tfc^r 1 = (I + T fe ,o + Ry 1 (I + fiF)- 1 

where /i — /i(A) is defined in $2.15 ), 



(3.2) sup ||iZ|| < C(Ko) |A| 

\k\<K 

for a positive constant C depending on go o,nd w, and 

F = S k P(I + T kt0 +R)- x . 
The operators R and F commute with complex conjugation. 
Proof. Let 

R = S k (A qx -A qo )Q. 
Since Q commutes with A 9A and A it follows that 

R=(S k Q) (A qx -A qo )Q. 

From Remark |2.3[ we have 

S k Q = S Q + H k Q 

so that sup| fe |<# \\S k Q\\ H -i/2^i\_y H i/2rgi^ is finite. We then use the fact that 
\W\ ~ ^-</o \\h 1 / 2 (s 1 )^h- 1 / 2 (s 1 ) — C \M 



to conclude that (3.2 1 holds. Write 



T k , x = nS k P + S k Q (A qx -A qn ) + S k Q {A qo - A ) 
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and use the fact that X^o = SkQ (A 9o — Ao) (since A qo P = AqP = 0) to conclude 
that (3.1) holds. The statement about complex conjugation follows from the fact 
that the projections P and Q, the operators A qx , and the operator Sk have the 
same property. □ 

Remark 3.2. The fact that Tk.o = SkQ (A 9o — Ao) implies that sup\ k \ <Ko ||Tfc,o| is 

bounded, and hence that sup| fc | <Jfo (I + Tk.o) 1 is bounded since k h-> Tu,o is 

a bounded continuous operator-valued function for \k\ < Kq and (I + Jfc ) o) 1 is 
known to exist for all k by |18) . 



The effect of Lemma [3T| is to focus attention on the rank-one operator F. If we 
write 

Ftp = (<p,ip)x 

a short computation shows that 

(<p, ■ )x 



(I + fxF) l =I 



while 



l + /x(v,x) 
D(k, A) := det (I + pF) =l + fx (v>, x) 



In our case we have, with ipo defined in (2.1 ), 
(3.3) <P=(I + T, 



fe,0 



R* 



(3.4) X =[Hi(k-)+log(\k\)}ip . 
We can prove: 

Lemma 3.3. For any Kq > and all \k\ < Kq and X sufficiently small, the function 
D(k,X) is real-valued, smooth in X and k for k ^ 0, and radial in k. Moreover, 
D(k,X) obeys the small-\X\ asymptotics 

(3.5) D(k, A) = 1 + (j,(X) (2irh - log \k\) + O (|A| |fc|) . 



Proof. From ( 3.3 >-( |3.4[ ) we compute (ip, x) = F\ + F%, where 

(3.6) A) = (2nh - log (ip , (I + T k , a + R)' 1 <^ ) 

(3.7) F 2 (k,X)=((l + T£ + ir)- 1 <p ,H l (k>)<po) 



In Q, h = iJi(O) = -7/2tt, and, in @, = iJi(z) - F x (0). Since 

(J + T fe .o + i?) _1 (^o = <A), 

we have 

F 1 (fc,A)=/x(A) (27rft-log|fc|) 
which is obviously radial in fc. To see that F2(k, A) is radial in k, we note that if 

Ck = Tk,a + R 

and Uj, is the unitary operator on H 1 / 2 {S 1 ) given by (U<f,f) (e 10 ) — /(e^ e+ ^), we 
have 

U,f,Ck = C^kU^ 

and 
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The second identity implies that H 1 {kz)^o = 'H-k'fio satisfies U<f, {H\{k o) ip (o)) = 
Hi(e l< ^k o)<po- Using these identities and the fact that U^ipa = ¥>o we easily deduce 
that F 2 (e^k, A) = F 2 (k, A) so that F 2 is radial in k. Finally, 



|F 2 (fe,A)| <C 



I + T, 



k,0 



R* 



sup \H(k -)\<C\k\, 

\Z\<1 



which gives the required asymptotics since |/i(A)| < C |A| for small |A|. The reality 
of D (k, A) follows from the fact that <po and H(k ■ ) are real-valued functions 
and that the operators T kj o and R commute with complex conjugation. Smooth 
dependence on the parameters k and A is a consequence of the smooth dependence 
of Tfe x on these parameters. □ 

Remark 3.4. It follows from the lemma that for |fc| < Kq and A sufficiently close 
to zero and k ^ 0, 

|V fc D(fe,A)| = J^i(l + 0(A)) 

is nonzero, so that the zero set 

Zx = {k e M 2 : k ^ 0, D(k, A) = 0} 

is locally a smooth curve. By radial symmetry, Zx is actually a circle of radius rx 
depending on A. We easily deduce from (3.5) and Lemma 2.4 that 



(3 - 8) r A ~«p rv . , 2V(Q)A 

Corollary 3.5. For A > small, the exceptional set £ is empty. 



Proof. We see from (3.8) that for fixed Kq and small positive A, rx > Ko and 
thus there are no singularities for \k\ < Kq. From Theorem 2.1 of [23], if |fc| > 
2 <7a||l~ 



C\\(l + |z| 2 ) 1//2 5A||i~ then the modified CGO solutions e 1/kz ip(z,k), exist, are 



unique, and are locally integrable. The function qx has support in \z\ < 1, therefore 
choosing Kq = sup 0<A<1 C||2<7a||l°° and using (1.3), there are no singularities for 
/c| > K . □ 

The purpose of the following lemma is to show that, even when there are no 
exceptional points of t\(k), this does not violate Nachman's result [TH Theorem 3] 
because t\(k) does not satisfy the small k decay requirement. 

Lemma 3.6. If (I + Tk t \) exists and is bounded for small k ^ and ^t(A) ^ 
then 

t ^ = -fo|M + 0(fc) 

Proof. We expand e tkz and e ikz for small k as 1 + C(|fc|) and then calculate (/ + 
T/ca)" 1 explicitly. 



tx(k) 



Si 



? (A qx -A )(I + T k ,x) 



— 1 Akz 



dS 



= / (l + 0(fc))(A 9A -A )(I + T k ,xr\l + O(k))dS 



(3.9) 



5'i 



(A 9A -A )(I + T k ,x)- 1 ldS + O(k) 
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Expanding (I + T/-)l using ( |2.8[ ) gives us 

(J + T fc Jl= [/ + (5b + «fc-(log|fc|)P)(A, A -A )]l 
= 1 + /*(A)(27r/i — log 



and so (I + Tk^^l = l/[l + fx{X)(2nh - log |fc|)]. Applying this to gives 

2?r^(A) 



tx(k) 



l + /i(A)(27r/i-log|A;|) 



2tt 



□ 



Fo r small A ^ 0, we see the assumptions of the le mma are satisfied. Remark 
(3.2) gives us that (I + Tk,\)~ 1 is bounded, and lemma (2.4) gives us that /x(A) ^ 0. 
Therefore we have that 



|tA(fc)| 



-2tt 



O(k) 



>c\k\ t 



log |*| 

for all c, e > and fc small, which shows that t\(k) does not decay fast enough for 
Nachman's Theorem 3 to apply. 
Next, we show: 

Lemma 3.7. There is a Xq > so that, for all A with < —A < Ao and any 

k c e Z\ 

lim |tx(*)| = oo. 

Proof. From the formula 

^,*) = (/ + T M r 1 (e ife W)(z) 

we compute 

(3.10) i>(z, k) = [(I + T fc , + R)- 1 (I + ^F)- 1 (e ifc «»)] (z) 



(7 + T k , a + R) [I - 



MA) 
D(k,X) 



(J + T M + #) X (z) + i?(z,fc) 



D(fc,A) 

where i?(z, fc) is regular in z, k. If A < and ko e Za, the zero set of D(k, A), we 
have 



lim 

fe— ¥k a 



D(k,X) 

Formula (1.12) and the fact that /x(A) is nonzero for small nonzero A imply that, 
to show 

lim |tx(*)| = oo, 



it suffices to show that 
(3.11) 

and 



lim inf 

fc->fe c 



> 



(3.12) 



lim inf 



Jikz 



(A g -A )(/ + T M + J R)- 1 x (z)dS(z) 



> 
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where ip and \ are given respectively by (3.3) and (3.4). Now 

<P= (I + Tkfi+R)' 1 

, -l 



Writing 



using the fact that 



■ T fe*o) 1 <Po = <^o 

T fc,oV?0 = ( A 9o - 
= ( A 9o " 



pn 

Pa+O (A) 



A ) QSkPo 
A ) Q(H(k ■ )) 



(see (2.8)) satisfies 



F fe > || <c|fc|, 

and using the fact that |fc c | ~ Ce c / A , we conclude that, for A < small, <p — 
ifo + O (A), that e lh " z ~ 1 + O (e c / A ) (c is a positive constant), and hence that 



(3.11) holds. 



To prove (3.12), first note that 

X =(H(k-)+\og\k\)p 
■0(A). 



1 



MA) 



Wo(e c/A ) 



since H is smooth and ~ e c / A . Next, note that (I + Tk.o + R) 1 (fio = ipo so 
that finally 



(A, - A ) (/ + T M + i?) 1 xl (z) d5(z) 



> 



O e 



,c/A 



which shows that (3.12) holds. This proves the lemma. 



□ 



Proof of Theorem 1.2 We have already shown that, for all sufficiently small neg- 
ative A, the scattering transform is singular on a circle of radius r\ with the 



asymptotic behavior (3.8), and for small positive A the set £ is empty by Corollary 



3.5 The behavior of t^(fc) near k — is given by Lemma |3.6| It remains to show 
that tA is smooth elsewhere. This follows from the fact that D(k, A) is smooth and 
nonzero away from the singular circle, the formula (3.10), and the explicit formula 
for t A . □ 

4. Computational methods 
4.1. Evaluating eigenvalues of DN maps. Given a radial potential q, we wish 



to compute numerically the eigenvalues [i n (q) defined in (2.13). The straightforward 



approach would be this: use the finite element method (FEM) to solve the Dirichlet 
problem (1.6) with f — p n . Then evaluate A q ip directly from (1.7) by numerical 



differentiation of the FEM solution. We will actually compute fio{l) in this way, 
but for fi n (q) with n ^ we can avoid the instability of numerical differentiation 
as explained below. 
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Consider the Neumann problem 
(4.1) (-A + q)u = in fi, 



du , 
dv 



\dQ = .9, 



where the mean value of g e H~ l / 2 (dVL) is zero. Define the Ncumann-to-Dirichlet 
map by 

K q g := u\an, 

where the condition f^ulon ds = 0. The functions ip n are eigenfunctions for the 
ND map: U q tp n = v n (q)ip n . 

Taking f — ip n with n =/= in ( |1.6| results in 

(4.2) (-A + q)u n = in O, u n \en = <p n , 

and we know that 



du„ 
dv 



\on = Aq^n = Hn{q)tPr. 



Consider the Neumann problem 

(4.3) (-A + q)u n = in Q, 



dv 



\an = Hn{q)<fr, 



Now since the solution u n is the same in (4.2 ) and (4.3 1 we see that lZ q (/j, n (q)tp n ) — 
ip n . If fJ, n (q) 7^ 0, by linearity we get the following connection between the eigen- 
values of the DN and ND maps: 

v n (q)¥n = K q (f n = — 

Vn(q) 



(4.4) 



Equation (4.4) provides us with a stable way to compute fi n (q) using the finite 
element method for the solution of (4.1), since there is no numerical differentiation 
involved in the evaluation of v n (q). 

4.2. Solution of the boundary integral equation. We explain how to solve 



equation (1.9) approximately by numerical computation. We follow the method 



described in [12] . The trick is to write the integral equation approximately as a 
matrix equation on the truncated Fourier series domain. 

Choose N > 0. We represent a function / e H s (dfl) approximately by the 
truncated Fourier series vector 

f(-N) 

n-N+i) 



f 



/(0) 



W-i) 

where the Fourier coefficients are defined for —N<n<Nby 

r 2n 



f(n) ■= (/, ¥>„) = 



1 

2tt 



d6. 
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By standard Fourier series theory we get for well-behaved / 

N 
n=-N 

Sk (Aq — Ao) using a matrix acting 
\n\(p n , so the (27V + 1) x (27V + 1) matrix 



Our goal is to approximate the operator T 
on the truncated Fourier basis. 

We know analytically that Ao<p n = |n|y n , so the (27V 
representing the operator Ao is 

(4.5) L := diag[iV, TV - 1 . . . , 2, 1, 0, 1, 2, . . . , TV - 1, TV]. 

Likewise, the DN map A q can be represented by a diagonal matrix containing the 
eigenvalues fi n (q) defined in (2.131: 

(4.6) 



j q := diag[fi N {q),fi N - 1 (q),. . ., /ii (q),fXo(q), Hi{q) 

Note that in (4.6) we made use of the symmetry (2.14|. 
By Lemma 2.2 we can write S k 



,H N - 1 (q),fji N (q)}. 



So + Wk — (log \k\) P. In our case of fl being 



the unit disc, the standard single layer operator So has the matrix 

1 1 . . 1 1 1 , 



So 



1 



i.i.o.i.i, 



TV — 1 TV 



Furthermore, the projection operator P defined in (2.3) can be represented by 

P = diag[0,...,0,l,0,...,0]. 

It remains to find a matrix Hj, for the operator T-L k . We define the elements of 
H fe = [H fe (m,n)] by 

(4.7) H fe (m,n) := (Hk<Pn,<pv 



1 

2^ 



271 



(H k e l 



dO. 



Here m € {—N, . . . , N} is the row index and n € {— N, . . . , N} is the column index. 
The function T-Lke 1 can be evaluated numerically by applying a quadrature rule 
to the integral in (2.9). For this we need to be able to compute point values of 
H k (z). By ([2J]) and (|2.6| and we can write 



(4.8) 



H k (z)=H 1 (kz) 



e^ z gi (kz)-G (kz) 



2ir 3 v ' v ' 2tt 

Now the evaluation of H k (z) is reduced to computing Faddeev's fundamental solu- 
tion gi(z), since everything else is explicit in the right hand side of (4.8). Following 
[H (3.10)], that can be done simply using formula 

(4.9) 9l (z) = -U-"Re(Ei(iz)), 

where Ei stands for the exponential-integral special function whose implementation 
is readily available in mathematical software packages. As explained in |21j . one can 
avoid evaluating the functions gi(z) and Gq(z) in (4.8) near the singularity at z = 
by calculating the harmonic function H k (z) first on a circle \z\ = R enclosing the 
evaluation domain, and then using the classical Poisson kernel to calculate H k (z) 
for \z\ < R. 

Approximate solution of (1.9) is now given in the frequency domain by 
^ = [/ + S fc (L 9 -L )]- 1 (e' S ^ 1 ) 



(4.10) 



= [1+ (S + H fc - (log \k\)P)(L q - LoTHe^lon), 
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where we used the decomposition (2.8), and e ikz \gn stands for the Fourier expansion 
of e lkz , calculated as follows. Write z = e l9 and compute as in [TU1 Section 2] to 
get 

(ik) n ^ n 

^— f— , n > 0, 



e lkz = a n(k)e me with a n (k) = 

n — — oo 

The vector e lkz \gQ thus takes the explicit form 



0, n < 0. 



e lkz \dn 








1 

?'/.' 

'.2 



Recall now the infinite-precision formula 



k 2 /2 



{ik) N /N\ 



t(k)= / e lkz (A q -A M-,k)ds. 
Jon 



Once the Fourier coefficient vector V'ln is solved from (4.10), set 



and define a function g : dQ 



9 = (L g - Lq^Iq 

C using the truncated Fourier series inversion: 



N 

E 

i=-iV 



g{ri)ip n {6). 



Then the scattering transform can be computed approximately using the formula 

r2n 



(4.11) 



t(fc) 



a ik exp( — i& 



> g (6) do. 



o 



The approximation in ( 4.11[ ) is most accurate for k near zero. 

Of course, the accuracy in (4.11) can be improved also by increasing N, but 



there is a limit to that. With reasonable computational resources it is possible 
to compute accurately enough up to \n\ < N — 16, but computation for 

\n\ > 16 quickly becomes an extremely hard problem. The reason is that the 
difference \/J, n (q) — \n\ \ becomes exponentially small as \n\ grows. We remark that 
one can achieve higher accuracy by computing the difference |/^ T1 (g) — |n|| directly 
analogously to the approach in [8], but in this work there was no need for that. 

5. Computational results 
We will study two examples numerically: the simple case q\ = Xw in Section |5.2| 



and a more complicated case q\ = Qq + Xw in Section 5.3 Here go is a nontrivial 
conductivity-type potential. 




Figure 1. Test function w(z) = w(\z\) defined by formula (5.1) 
with Ri = 0.8 and R x = 0.9. 

5.1. Definition of a test function w. We define a radial Cg function w(z) — 
w(\z\) as follows. Take two radii < R\ < R 2 < 1 and define w(|z|) in three pieces: 

for < \z\ < 

(5.1) w(\z\)={ p(\z\) for R 1 <\z\<R 2 , 

for R 2 < \z\ < 1. 




The polynomial p in (5.1) is constructed as follows. Note that the polynomial 



p(t) = 1 — lOt + 15t 4 — 6t is smooth in the interval < t < 1 and satisfies 
p"(0) = p'(0) = = p'(l) = p"(l). Set for R 1 <t<R 2 

The test function defined above is twice continuously differentiable instead of in- 
finitely smooth as in the theoretical part above. However, the discrepancy is not 
essential in this numerical work. See Figure [T] for a plot of the test function w. 

5.2. First example: zero potential at A = 0. We set q\ = Xw with the radial 



test function defined by (5.1) with Ri = 0.8 and Ri = 0.9. 

Our aim is to compute the radial scattering transform t^(|A:|) for < |fc| < 3.5 
and for the parameter A ranging in a suitable interval. In practice we choose the 
following finite set of fc-values: 

(5.2) k = 0.01, 0.02, 0.03, . . . , 3.49, 3.50. 

Note that the fc-grid is bounded away from zero by a significant gap of size 10~ 2 . 
Furthermore, we consider the following choices of parameter A: 

(5.3) A = -35.00, -34.95, -34.90, . . . , 34.90, 34.95, 35.00. 



We start the numerical computations by constructing the matrices (4.7) for each 



fc-value listed in (5.2 1. We take N = 12, so each Hj, has size 25x25. These matrices 
need to be computed only once for a given fc, so we can reuse the matrices in our 
second example below. 
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A A 

Figure 2. Eigenvalues //(A) = ^o(Qx) corresponding to the first 
example potential. Left: plot with full parameter range — 35 < 
A < 35. Right: detail of the left plot. Note that fi(0) = and 
£t'(0) > as predicted by Lemma [2~4{ Also, note that Dirichlet 
eigenvalues of the potential q\ in the unit disc cause singularities 
in //(A). 



Next we use the methods described in Section |4.1| to compute the eigenvalues 
of the dn map corresponding to each potential q\ with A ranging as in (5.3 1. We 
construct a finite element mesh for the unit disc with 131585 nodes and 262144 
triangles. We use Matlab's PDE toolbox to solve the Neumann problem (4.1 1 with 
g = if n for 1, . . . ,N, and get accurate approximations to the eigenvalues v\, . . . , 
of the nd map. By (2.14) and (4.4) we see that we know all eigenvalues n n {q\) of 
the DN maps except for /j,o(q\). We use FEM to solve Dirichlet problems of the 
form ( |1.6[ ) where q = q\ and / = 1; this way we get good approximations to ^o(qx). 
See Figure|2]for plots of the eigenvalue Ho(qx) as function of the parameter A. 

Now we have all the ingredients for solving the matrix equation (4.10) and eval- 
uating the scattering transform by (4.11). We used Matlab, a parallelization mid- 
dleware solution provided by Techila Ltd, and two hardware solutions: the "Ukko" 
cluster computer of the Computer Science Department of University of Helsinki, 
and Microsoft Azure cloud computing services. 

See Figure [3] for a plot of the scattering transform profiles as a grayscale plot, 
and Figure [4] for some selected profiles as conventional plots. Further, see Figure 
[5] for a comparison of numerical results and the asymptotic formula ( |1.5[ ) for the 
radius of the exceptional circle. 
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A 

Figure 3. Scattering transform corresponding to the first ex- 
ample. The horizontal axis is the parameter A in the definition 
q\(z) — Xw(z) of the potential. The vertical axis is |fc|. There are 
curves along which a singular jump "from — oo to +00" appears. 
The k values at those curves are exceptional points. See Figure [4] 
for further illustration of the singularities. 




1*1 

Figure 4. Profiles of scattering transforms for various A. This 
picture corresponds to Example 1. 
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A A 
Figure 5. Comparison of numerical results and the asymptotic 



formula (1.5| for the radius of the exceptional circle. This plot 
is for Example 1. Left: detail from Figure [3] with parameters 
ranging in the rectangle -4 < A < and 0.001 < |fc| < 0.4. 



Right: the asymptotic function r(A) given by Theorem 1.2 For 
ease of comparison, we also show in the background the pixel image 
from the left but with a lighter colormap. The asymptotic formula 
matches the computational result very closely in the interval — 2 < 
A < 0. 
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5.3. Second example: nontrivial conductivity-type potential at A = 0. 

Here we define q\ — tfo + Au>, where the test function w is defined by formula (5.1 ) 
and plotted in Figure [l] The potential §o corresponding to A — is defined by 
q = a-" 1 / 2 Act 1 / 2 . See Figure for plots of the conductivity a and the potential 
§o- Figure [7] shows the profile of the non-singular scattering transform of §o- 

We compute the eigenvalues of the dn map corresponding to each potential q\ 
similarly than in Example 1. See Figure[8]for plots of the eigenvalue fi(X) = Ho(q\) 
as function of the parameter A. 

We compute the scattering transform similarly to Example 1 above. See Figure 
[9] for a plot of the scattering transform profiles as a grayscale plot. Further, see 
Figure 10 for a comparison of numerical results and the asymptotic formula (1.5) 
for the radius of the exceptional circle. 
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Conductivity a(z) Profile of conductivity 




Conductivity-type potential go (2) Profile of potential 




?o(M) 



-1 -1 



Figure 6. Top row: mesh plot and profile plot of the rotationally 
symmetric conductivity <j(z) = a{\z\). Bottom row: mesh plot and 
profile plot of the resulting conductivity- type potential cfo(z) = 

?o(M). 
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Scattering transform to(fc) 



Profile of scattering transform 



0.5- 
0- 
-0.5- 
-1- 
-1.5, 

-2-5-1 
3 




0.5 


-0.5 
-1 
-1.5 

-2 
-2.5 



to(|fc|) 



Figure 7. Left: mesh plot of the rotationally symmetric scatter- 
ing transform to(fc) = to(|fc|) corresponding to the initial potential 
§o shown in Figure [6] Right: profile plot of t (|fc|). 




A A 

FIGURE 8. Eigenvalues fi(X) = Ho(q\) corresponding to the second 
example potential. Left: plot with full parameter range —35 < A < 
35. Right: detail of the left plot. Note that ^(0) = and //(0) > 



as predicted by Lemma 2.4 Also, note that Dirichlet eigenvalues 
of the potential q\ in the unit disc cause singularities in /i(A). 
Compare to Figure [2] 
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35 



25 



1*1 



1.5 



0.01 













-35 -30 -25 -20 -15 -10 -5 5 10 15 20 25 30 35 



A 



Figure 9. Scattering transform corresponding to the second ex- 
ample. The horizontal axis is the parameter A in the definition 
Qx = Qo + Aw of the potential. The vertical axis is There are 
curves along which a singular jump "from — oo to +00" appears. 
The k values at those curves are exceptional points. Compare to 
Figure [3} 
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Figure 10. Comparison of numerical results and the asymptotic 
formula (1.5 1 for the radius of the exceptional circle. This plot 
is for Example 2. Left: detail from Figure [9] with parameters 
ranging in the rectangle —4 < A < and 0.001 < \k\ < 0.6. 
Right: the asymptotic function r(A) given by Theorem 1.2 For 
ease of comparison, we also show in the background the pixel image 
from the left but with a lighter colormap. The asymptotic formula 
matches the computational result very closely in the interval — 2 < 
A < 0. Compare to Figure [5] 
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6. Discussion 

We study zero-energy exceptional points of radial, compactly supported Schrodinger 
potentials in dimension two. Our work was inspired by preliminary numerical 
experiments showing the emergence of singularities in the scattering transforms of 
potentials constructed by subtracting a test function from a radial conductivity- 
type potential. 

We prove new results for radial, real-valued potentials of the form q\ = q + Xw, 
where q is of conductivity type, the test function w is non-negative, and A el. It 
turns out that for small positive A there are no exceptional points and for small 
negative A there is exactly one exceptional circle whose radius we can determine 
asymptotically. In our two computational examples the asymptotic formula is quite 
accurate in the rather substantial interval — 2 < A < 0, see Figures [5| and [TO] See 
Figure [TTJ for a diagram illustrating all currently known theoretical facts about 
two-dimensional exceptional points at zero energy. 

Our numerical computations raise some further theoretical questions as well. 
Figures [3] and [9] suggest that for negative A far away from zero there are several 
exceptional circles. Also, for large positive A some exceptional points appear, but 
it is unclear where in the (A, |fc|) plane they originate. Furthermore, there are some 
quite complicated features in the parameter interval —23 < A < —18.5 that are 



very different for our two examples. See Figure 12 for a zoom-in 



Appendix A. Radial symmetry of scattering transforms 



Let q\ be a potential of the form (1.4). If k G C \ is not an exceptional point 
of q\, then uniqueness of the CGO solution ip(z,k) shows that all k! £ C with 
|fc'| = |fc| are non-exceptional for q\ as well. Furthermore, we can argue as in |14[ 
Section 4.1] and find out that the scattering transform satisfies t\(k) = t a ( | -fc | ) and 
txjkj = t x (k). 

Also, we know from Remark |3.4| that exceptional points appear on a set of 
measure zero (circles centered at the origin). So, for illustrating t\(k) it is enough 
to display real-valued profiles of t^(|A;|) evaluated at |fc| > 0, as is done in Figures 
land[9l 

Appendix B. Conductivity- type potentials and criticality 

Let q(z) = q(\z\) be a radial potential of conductivity type in the sense of Definition 
|l~l~j Then we can write q — tp 1 (Aip) with some smooth and strictly positive 
function tp for which ^ — 1 is compactly supported. Set 

q\ = q + Aw 

with w <E C^°(M 2 ) a radial, nonnegative test function which is not identically zero. 

First of all, [TTl Theorem 3.1(iii)] implies that qo is critical. To see this, note 
that ffo(M) = VKM) is the unique solution of [TTl equation (3.3)] with j = and 
n = 2 satisfying the asymptotic condition g (\z\) — 1 + o(l) as \z\ —> oo. Then the 
integral in [T7J formula (3.5)] diverges. 

It follows from [T7J Theorem 2.4(i)] that the potential q\ is supercritical for all 
A < 0. In that case we know from [TJ [TBI [5] that there is no positive solution of 
(—A + q\) ij) — 0, so q\ is not of conductivity type. 

It follows from [TTl Theorem 2.5(i)] that the potential q\ is subcritical for all 
A > 0. Furthermore, by [TTJ Theorem 5.6(i)] the unique positive solution of 



No exceptional points 



9j 




\ 



-An An 



Figure 11. Conceptual diagram illustrating the current under- 
standing of exceptional points in dimension two. The underlying 
grayscale image shows the scattering transform of Figure [3] corre- 
sponding to Example 1. The potential go is of conductivity type, 
and therefore by Nachman |18j there are no exceptional points for 
A = 0. Theorem 1.2 above implies that there is a An. > such 



that (i) for < A < A the potential q\ has no nonzero exceptional 
points, and (ii) for — An < A < the complex numbers satisfying 
|fc| = r(A) are the only nonzero exceptional points for q\. Fur- 
thermore, the Neumann series argument in the proof of [181 Thm 
1.1] shows that there exists aradius R(X), depending on ||<?a|!lp(r 2 ); 
such that q\ does not have exceptional points satisfying |fe| > R(X). 
Note that the number An and the radius i?(A) shown above are not 
in scale. 



(— A + q\)ip = has asymptotics c log | z | + O(l) with a positive constant c as 
\z\ — > oo. Thus q\ is not of conductivity type because that would require a finite 
and constant limit limi z |_>oo 
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Example 1 Example 2 




-23 -22 -21 -20 -19 -23 -22 -21 -20 -19 



A A 

Figure 12. Comparison of the scattering transforms of the two 
example potentials in the parameter domain —23 < A < —18.5 and 
< |fc| < 1.5. 

Appendix C. Spectral Properties of Conductivity-Type Potentials 

The purpose of this appendix is to analyze the spectrum of the self-adjoint 
operator H — —A + q when q is a potential of conductivity type. 

Proposition C.l. Suppose that qo G Cg°(R 2 ) is a potential of conductivity type, 
and let H = —A + qg. Then H has no L 2 -eigenvalues. 

Proof. By standard arguments (see, for example, Theorem XIII. 56 in [20), the 
equation Hip = A?/> has no L 2 solutions for any A > 0, while the positivity of the 
quadratic form associated to H shows that the there can be no eigenvalues with 
A < 0. It remains to show that, also, there are no solutions of Hip = that vanish at 
infinity, hence no L 2 -eigenvalues at zero energy. This is an immediate consequence 
of arguments in Nachman |18j but we reproduce them for the reader's convenience. 

Suppose that h e H 1 ^ 2 ) is a weak solution of Lh = Without loss we may 
assume that h is real-valued. By elliptic regularity, h is a bounded, C°° function. 
Let 

v = hdipo — Tpodh. 

Note that, as ip Q £ L°°, d%p a € C§°, and h e H 1 , it follows that v e L 2 . A 
straightforward computation using the facts that Addtpo — qipo and Addh = qh 
shows dv = av — av where a = dipo/ipo. Note that a 6 Cg°. By a standard vanishing 
theorem for generalized analytic functions (see for example [28: , we conclude that 
v = 0, and hence that d (Ji/tpo) = 0, i,.e., h/tpo is antiholomorphic. Since h vanishes 
at infinity while -00 is bounded below, we conclude that h = cipo for a constant c, 
and hence, h = since ipo{z) — >• 1 as \z\ — > oo. □ 

Remark C.2. The argument above shows that a conductivity- type potential is rep- 
resented by a unique normalized positive solution tpQ. 

Next, we consider spectral properties of — A + q on a bounded domain containing 
the support of q. Recall that Hq(H) is the completion of C^(£l) in the H^-norm. 
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Proposition C.3. Let Q. be a bounded domain in M 2 with smooth boundary, so 
chosen that suppgo * s strictly contained in fi. Denote by Hq the operator — A + q 
with Dirichlet conditions on dfl. Then is not a Dirichlet eigenvalue of Hq ■ 

Proof. Suppose that ip £ Hq(SY) satisfies Hntp = 0. Extend ip to a function \ m 
i? 1 (]R 2 ) by setting x( z ) = f° r z £ K 2 \^- The distribution gradient of x is given 
by (Vx) ( z ) — {^i } ){ z ) f° r z £ an d (Sx){ z ) = otherwise. Letting q be the 
quadratic form of —A + qo on H 1 (M 2 ), we have 

q(x,x) = o 

so that, for any ip £ C§°(M 2 ), the function 

F(t) = q(x + t<P,X + t<p) 
has an absolute minimum at t = 0. Since F'(0) — we recover Req(tp,x) = 
for any such ip. It follows that q(<p,x) = f° r & U <P G Co°(K 2 ), hence x is a weak 
solution of Hx = which vanishes identically outside f2. We can now use unique 
continuation arguments (see, for example, |20j . Theorem XIII. 57) now show that 
X = 0, hence ip = 0. □ 
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